Coverage Probability of Wald Interval for Binomial 

Parameters * 



Xinjia Chen 
Submitted in April, 2008 

Abstract 

In this paper, we develop an exact method for computing the minimum coverage probability 
of Wald interval for estimation of binomial parameters. Similar approach can be used for other 
type of confidence intervals. 

1 Wald Interval 

Let X be a Bernoulli random variable with distribution Pr{X = 1} = 1 — Pr{X = 0} = p E (0, 1). 
Let Xi, • ■ • , X n be i.i.d. random samples of X. Let K = Ya=i The widely-used Wald interval 
is [L, U] with lower confidence limit 



L 




and upper confidence limit 



u = - + z 5/2 \^ — m 

n V n 

x 2 

where Z$/2 1S the critical value such that ^ -j=e~~ = |. It has been discovered by Brown et 
al. p] that the coverage probability of Wald interval is surprisingly poor. 



2 Coverage Probability 

The coverage probability of Wald interval for binomial parameters has been investigated by [TJ and 
other researchers by Edgeworth expansion method and numerical methods based on discretizing 
the binomial parameter. Here, we have obtained expression of the minimum coverage probability 
of Wald interval, which requires only finite many evaluations of coverage probability. 
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Theorem 1 Define 

_ 2 P + e-^e^ + Ae P {i-p) _ 2 P + 6 + yg + A6 P (i- P ) 

{P) 2(1 + 6) ' [P> 2(1 + 6) 

z 2 

with 6 = Define 

u I ki\ _ k) u I k(\ _ k\ 

J?(k) = -- Z 5/2 \ si < W(k) = - + Z 5/2 \ ^ *L 

n ' V n n V n 

for k = 0, 1, • • • ,n. Define 

Q(k) = Pr{pr-(jSf(jb))l <K<k-l\ JSf(fc)}, C[(k) = Pr{LT - (j2f(fc))J + 1 < K < k - 1 | Sf(k)} 

for k G {0, 1, • • • ,n} suc/i £/ia£ < S£(k) < 1. Define 

C u (k) = Pr{fc + 1 < x < LT+(^(fc))J | ^r(fc)}, C4(fc) = Pr{fc + 1 < A' < \T+(W(k))] - 1 | 

fork G {0, 1, • • • ,n} suc/t i/iai < <2r(fc) < 1. 

Su PP ose 6 < 3. Then, the following statements hold true: 
(I): inf pg (o,i) P r {-^ < P < U \ p} equals to the minimum of 

{Q(k) : < k < n; < Jzf(fc) < 1} U {C u (k) : < fc < n; < ^(fc) < 1}. 

(TJj: min pg (o,i) Pr{^ < P < f I p} equals to the minimum of 

{C[(k) : < k < n; < jSf(jfc) < 1} U {C^(fc) : < fc < n; < <%(k) < 1}. 



The proof of Theorem 1 is given in the next section. 

3 Proof of Theorem 1 

We need some preliminary results. 

z 2 

Lemma 1 For n > — ^p-, both the lower and u PP er confidence limits of Wald interval are mono- 
tonically increasing with res P ect to k. 

Proof. For simplicity of notation, let 

_ k 
n 

Then, the upper confidence limit can be written as 

h(z) = z + Z S /2\/z(l - z)/n. 
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Similarly, the lower confidence limit can be written as 

g(z) = z - Z s/2 y/z(l - z)/n. 

To show that the upper confidence limit is monotonically increasing with respect to k, it suffices 
to show that 

OZ 

if < h{z) < 1. Since 

dh(z) \/9 l-2z 



1 + 



dz ~ 2 yfi(T=Z)' 

which is clearly positive for < z < |, it remains to show 
or equivalently, 

z{l-z)>-(2z-l)\ Vzef^l 



4' ' V 2 

Note that h(z) < 1 for < z < j^, and h(z) > 1 for 1 > z > j^. 

If < ~, i.e., 9 > 1, then we are done, since > for all < 2 < < ~. Otherwise, 
if 9 < 1 , it suffices to show 

a 

w(z) = z(l - z) - -(2z - l) 2 > 

for | < z < j^tq. Since z(l — z) is decreasing and \{2z — l) 2 is increasing for 1 > z > |, is 
decreasing for ^ < z < Therefore, it suffices to have 



w ( IT 77 1 > °' 



i.e., 



1 + V 1 + / 4 V 1 + 

i.e., 

1 ^ 2 



2 

1 ) > o, 



i.e., 



(1 + 0) 2 4V1 + 



4 > (1 



>0, 



which is guaranteed since < 9 < 3. This shows that the upper confidence limit is monotonically 
increasing with respect to k. Observing that g(z) = l — h(l — z), we have that the lower confidence 
limit is also monotonically increasing with respect to k. 

□ 
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Now we consider the minimum coverage probability. By solving equation 




with respect to k, we can show that 



Pr{L <p<U 




= Pr{k < K < T+0) 


W(k)}, 


< ^(jfe) < 1 


Pr{L <p<U 


J2f(fc)} 


= Pr{T"(p) <K <k 


J2f(fc)}, 


< Jgf (fc) < 1 


Pr{L < p < U 


W(k)} 


= Pv{k < K < T + (p) 


3r(fc)}, 


< ^(jfe) < 1 


Pi{L < p <U 


<?(k)} 


= Pr{T" (p) <K <k 


^(fc)}, 


< Jgf (jfe) < 1 



Since both the lower and upper confidence limits of Wald interval are monotone as asserted by 
Lemma 1, the proof of Theorem 1 can be completed by making use of the above results and 
applying the theory of coverage probability of random intervals established by Chen in [2]. 
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